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We study graphene-based electron systems with long-range Coulomb interaction by performing 
an analytic continuation in the number of dimensions. We characterize in this way the crossover 
between the marginal Fermi liquid behavior of a graphite layer and the Luttinger liquid behavior at 
D = 1. The former persists for any dimension above D = 1. However, the proximity to the D = 1 
fixed-point strongly influences the phenomenology of quasi-onedimensional systems, giving rise to 
an effective power-law behavior of observables like the density of states. This applies to nanotubes 
of large radius, for which we predict a lower bound of the corresponding exponent that turns out to 
be very close to the value measured in multi-walled nanotubes. 
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Recently there has been much interest in the search 
of unconventional electron behavior deviating from the 
Fermi liquid picture JjJ. Besides this, the other paradigm 
that is well-established on theoretical grounds is the Lut- 
tinger liquid behavior of one-dimensional (ID) electron 
systems There have been suggestions that this be- 

havior could be extended to two-dimensional (2D) sys- 
tems, in the hope that it may explain some of the fea- 
tures of the copper-oxide materials Q. However, at least 
for the Luttinger model, the analytic continuation in the 
number D of dimensions has shown that the Luttinger 
liquid behavior is lost as soon as one departs from D = 1 

Several authors have also analyzed the possibility that 
singular interactions could lead to the breakdown of 
the Fermi liquid picture Q. With regard to real low- 
dimensional systems, such as carbon nanotubes, the main 
electron interaction comes actually from the long-range 
Coulomb potential V(|r|) ~ 1/ |r| . This is also the case of 
the 2D layers in graphite, which have a vanishing density 
of states at the Fermi level. Quite remarkably, a quasi- 
particle decay rate linear in energy has been measured 
experimentally in graphite 0, pointing at the marginal 
Fermi liquid behavior in such 2D layers. Due to the singu- 
lar Coulomb interaction, the imaginary part of the elec- 
tron self-energy in the 2D system behaves at weak g cou- 
pling like g 2 u> M . It is crucial, though, the fact that the 
effective coupling scales at low energy as g ~ l/log(w). 
This prevents the logarithmic suppression of the quasi- 
particle weight, which gets corrected by terms of order 
5 2 logH~l/logH @. 

In this letter we investigate whether the long-range 
Coulomb interaction may lead to the breakdown of the 
Fermi liquid behavior at any dimension between D = 1 
and 2. The issue is significant for the purpose of compar- 
ing with recent experimental observations of power-law 
behavior of the tunneling conductance in multi-walled 
nanotubes |ll[] . These are systems whose description lies 
between that of a pure ID system and the 2D graphite 



layer. It turns out, for instance, that the critical exponent 
measured for tunneling into the bulk of the multi-walled 
nanotubes isa~ 0.3. This value is close to the exponent 
found for the single- walled nanotubes p2 13 1. However, 
it is much larger than expected by taking into account 
the reduction due to screening (~ l/\/~N) in a wire with 
a large number TV of subbands, what points towards sen- 
sible effects of the long-range Coulomb interaction in the 
system. 

We develop the analytic continuation in the number 
of dimensions having in mind the low-energy modes of 
metallic nanotubes, which have linear branches crossing 
at the Fermi level. From this picture, we build at general 
dimension D a manifold of linear branches in momentum 
space crossing at a given Fermi point. We consider the 
hamiltonian 
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where the <ii matrices are defined formally by {at, <jj} — 
2Sij. Here p(p) are density operators made of the elec- 
tron modes V&(p), and c(Z?)/|p|' D_1 corresponds to the 
Fourier transform of the Coulomb potential in dimen- 
sion D. Its usual logarithmic dependence on |p| at 
D = 1 is obtained by taking the ID limit with c(D) = 
r((D-l)/2)/(2^?- D . 

The dispersion relation e(p) = ±|p| is that of Dirac 
fcrmions, with a vanishing density of states at the Fermi 
level above D = 1. This ensures that the Coulomb inter- 
action remains unscreened in the analytic continuation. 
At D — 2 we recover the low-energy description of the 
electronic properties of a graphite layer, dominated by 
the presence of isolated Fermi points with conical disper- 
sion relation at the corners of the Brillouin zone . 

In the above picture, we are neglecting interactions 
that mix the two inequivalent Fermi points common to 
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the low-energy spectra of graphite layers and metallic 
nanotubes. In the latter, such interactions have been 
considered in Refs. |L5| and p6| , with the result that 
they have smaller relative strength (~ 0.1/-ZV, in terms of 
the number TV of subbands) and remain small down to 
extremely low energies. More recently, the question has 
been addressed regarding the interactions in the graphite 
layer, and it also turns out that phases with broken sym- 
metry cannot be realized, unless the system is doped 
about half-filling (TJ or it is in a strong coupling regime 

We will accomplish a self-consistent solution of the 
model by looking for fixed-points of the renormalization 
group transformations implemented by the reduction of 
the cutoff A jl9| . As usual, the integration of high-energy 
modes at that scale leads to the cutoff dependence of the 
parameters in the low-energy effective theory. We will see 
that the Fermi velocity vf grows in general as the cutoff 
is reduced towards the Fermi point. On the other hand, 
the electron charge e stays constant as A — > 0. This 
comes from the fact that the polarizability II does not 
show any singular dependence on the high-energy cutoff 
A for D < 3. The polarizability is then given by 
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where b(D) 
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The dependence of vf on the cutoff A implies an in- 
complete cancellation between self-energy and vertex cor- 
rections to the polarizability. The dressed polarizability 
depends therefore on the effective Fermi velocity t>^(A). 
The renormalized value of vp is determined by fixing 
it self-consistently to the value obtained in the electron 
propagator G corrected by the self-energy contribution 
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The fixed-points of the renormalization group in the 
limit A — > determine the universality class to which 
the model belongs. At D — 2, we are bound to obtain 
the low-energy fixed-point at vanishing coupling of the 
model of Dirac fermions with Coulomb interaction |[o[ . 
On the other hand, at D = 1 there has to be presumably 
a fixed-point corresponding to Luttinger liquid behav- 
ior. We note, however, that no solution of the model 
has been obtained yet without carrying dependence on 
the transverse scale needed to define the ID logarithmic 
potential. Our dimensional regularization overcomes the 
problem of introducing such external parameter, which 
prevents a proper scaling behavior of the model P0|] . 

At general D, the self-energy (||) shows a logarithmic 
dependence on the cutoff at small frequency u>k and small 
momentum k. This is the signature of the renormaliza- 
tion of the electron field scale and the Fermi velocity. In 



the low-energy theory with high-energy modes integrated 
out, the electron propagator becomes 
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where g = b{D)c(D)e 2 /v F , f(D) 

h n (D) — F p"^ 3°d) 2 /2+i) ■ ^ ne q uan tity Z 1 / 2 represents 
the scale of the bare electron field compared to that of 
the renormalized electron field for which G is computed. 

The renormalized propagator G must be cutoff- 
independent, as it leads to observable quantities in the 
quantum theory. This condition is enforced by fixing the 
dependence of the effective parameters Z and vp on A as 
more states are integrated out from high-energy shells. 
We get the differential renormalization group equations 
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At D = 2, the right-hand-side of these equations can 
be summed up to the functions that have been found 
previously in the renormalization of the graphite layer 
ETJ . Furthermore, they also provide meaningful expres- 
sions in the ID limit. At D = 1, the right-hand-side of 
Eq. (H) vanishes identically as a function of the variable 
g. Therefore, the ID model has formally a line of fixed- 
points, as it happens in the case of a short-range inter- 
action. The scaling of the electron wavefunction can be 
read from the right-hand-side of Eq. (^|) , which becomes 
(2 + .g)/(2 v / TTg) - 1 at D = 1. This coincides with the 
anomalous dimension that is found in the solution of the 
Luttinger model, what provides an independent check of 
the renormalization group approach to the ID system. 

We have therefore a model that interpolates between 
marginal Fermi liquid behavior, that is known to char- 
acterize the 2D model, and non-Fermi liquid behavior 
at D = 1. As the electron charge e is not renormal- 
ized for D < 3, the scaling of the effective coupling 
g = b(D)c(D)e 2 jvp is given after Eq. (0) by 
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The right-hand-side of Eq. (^) is a monotonous increas- 
ing function of g, for any dimension between 1 and 2, as 
observed in Fig. |l| . 
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FIG. 1. Plot of the P function at the right-hand-side of Eq. 
([?]) for different dimensions. From top to bottom, the curves 
correspond to the values D = 2.0, 1.7, 1.4, 1.3, 1.2 and 1.1 . 



We find that, away from D = 1, there is only one 
fixed-point of the renormalization group at g = 0. The 
scale dependence of the effective coupling e 2 /vp is dis- 
played for different values of D in Fig. 0, where the 
flow to the fixed-point is seen. Consequently, the scale 
Z of the wavefunction is not renormalized to zero in the 
low-energy limit, and the quasiparticle weight remains 
finite above D = 1. We conclude then that, even in a 
model that keeps the Coulomb interaction unscreened, 
the breakdown of the Fermi liquid behavior only takes 
place formally at D = 1. 




The subtlety concerning the long-range Coulomb in- 
teraction is that the function c(D) diverges in the limit 
D — > 1. This is actually what transforms the power-law 
dependence of the potential into a logarithmic depen- 
dence at D = 1. We observe that the ID limit and the 
low-energy limit A — * do not commute. If we stick to 
D = 1, we obtain a divergent coupling g for the Coulomb 
interaction as well as a divergent electron scaling dimen- 
sion. At any dimension slightly above D = 1, however, 
the fixed-point is at g = 0, with its corresponding van- 
ishing anomalous dimension. 

In order to understand whether the ID model has any 
stable fixed-point for finite values of e 2 /vf, one can study 
the model by performing an expansion in powers of g~ x ■ 
The Fermi velocity is renormalized by terms that are an- 
alytic near the point g = oo, and that lead to the scaling 
equation 

A-±MA) = v F f(D) U-D-l) T /2 ^\ (*) 
dA ' * J ^'y d J y((D + 1)/2) v ' 

up to terms of order 0(g^ 1 ). In the limit D — > 1, g — > oo, 
the right-hand-side of Eq. (H) vanishes identically. This 
confirms, on nonperturbative grounds, that the ID model 
with the Coulomb interaction has a line of fixed-points 
covering all values of e 2 /vf- 

In the vicinity of D = 1, the presence of such critical 
line becomes sensible, and a crossover takes place to a be- 
havior with a sharp reduction of the quasiparticle weight. 
This can be seen in the renormalization of the electron 
field scale Z, displayed in Fig. ||[ For values of D above 
w 1.2, we have a clear signature of quasiparticles in the 
value of Z at low energies. For lower values of D, the 
picture cannot be distinguished from that of a vanishing 
quasiparticle weight for all practical purposes. The dras- 
tic suppression of the electron field scale Z takes place 
over a variation of only two orders of magnitude in the 
energy scale. 




FIG. 2. From top to bottom, curves of the effec- 
tive coupling as a function of the energy scale for 
D = 1.1, 1.2, 1.3, 1.4, 1.7 and 2.0 . 



FIG. 3. Plot of the electron field scale Z for a value of the 
bare coupling e 2 /(n 2 VF) = 4.0. The different curves corre- 
spond, from top to bottom, to D — 2.0, 1.7, 1.4, 1.3, 1.2 and 
1.1 . 
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The above picture allows us to make contact with the 
experiments carried out in multi- walled nanotubes 
In the proximity of the D = 1 fixed-point, the density of 
states displays an effective power-law behavior, with an 
increasingly large exponent. Moving to the other side of 
the crossover, the density of states approaches the well- 
known behavior of the graphite layer, n(e) ~ |e|. In Fig. 
H we give the representation of the density of states 

n(e) ~ Zie)^ - 1 (9) 



corresponding exponent that turns out to be very close 
to the value measured experimentally. 
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for several dimensions approaching D = 1. 




FIG. 4. Plot of the logarithm of the density of states for 
a value of the bare coupling e 2 /(tt 2 vf) = 0.5. The solid line 
corresponds to D — 1.1 and the dashed lines correspond, from 
top to bottom, to D = 1.2, 1.3 and 1.4 . 

A value e 2 /(ir 2 VF) — 0.5 for the bare coupling is ap- 
propriate for typical multi-walled nanotubes, as it takes 
into account the reduction due to the interaction with the 
inner metallic cylinders pi] . We observe that the expo- 
nents of n(e) at different dimensions are always larger 
than a lower bound a ~ 0.26. This is in agreement 
with the values measured experimentally. Our analy- 
sis stresses the need of an appropriate description of the 
dimensional crossover between one and two dimensions, 
showing that the picture of a thick nanotube as an aggre- 
gate of ID channels does not allow to obtain the correct 
values of the critical exponents. 

To summarize, we have studied the renormalization of 
the Coulomb interaction in graphene-based structures. 
We have made a rigorous characterization of the differ- 
ent behaviors, as we have proceeded by identifying the 
fixed-points of the theory We have seen that the Fermi 
liquid behavior persists formally for any dimension above 
D = 1, as it also happens in the case of a short-range in- 
teraction ||. On the other hand, the proximity to the 
ID fixed-point influences strongly the phenomenology of 
real quasi-onedimensional systems, giving rise to an effec- 
tive power-law behavior of observables like the tunneling 
density of states. This is the case of the multi-walled 
nanotubes, for which we predict a lower bound for the 
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